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[^- I We consider strong and electromagnetic isospin violation in near-threshold pion-kaon 

scattering. At tree level, such effects are small for all physical channels. We work out 
the complete one-loop corrections to the process ir~K + — » tt°K°. They come out 
rather small. We also show that the corresponding radiative cross section is highly 



suppressed at threshold. 
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1 Introduction 

The quark masses allow to consider various approximations to Quantum Chromodynamics 
(QCD). For the c, b, and t quarks, the masses are so large that an expansion in inverse powers 
of these masses can be performed systematically. This leads to the so-called heavy-quark 
effective field theory. On the other hand, to a good first approximation, one can consider 
the light quarks u, d, and s as massless. In that limit, the QCD Lagrangian exhibits a 
chiral symmetry which is, however, spontaneously broken as witnessed by the appearance 
of eight almost massless pseudoscalar Goldstone bosons, the pions, the kaons, and the eta. 
These would-be Goldstone bosons acquire their masses from the explicit symmetry breaking 
due to the quark mass term. Spontaneous as well as explicit symmetry violation can be 
explored systematically by means of chiral perturbation theory (ChPT), an effective field 
theory formulated in terms of the asymptotically observable fields. At low energies, the 
generating functional of QCD is characterized by two energy scales. One is the pion (kaon) 
decay constant in the chiral limit, denoted by F . Its non-vanishing value (F ~ 88 MeV) 
is a necessary and sufficient condition for spontaneous symmetry breaking, much like the 
vacuum expectation value of the neutral Higgs boson signals the breaking of the electroweak 
symmetry. The second scale is given by the quark condensate, B = |(0|gg|0)|/.F 2 . In the 
standard scenario of chiral symmetry breaking, (0|gg|0) ~ (—225 MeV) 3 so that B ~ 
1.5 GeV^> F. This leads e.g. to a very precise prediction for the S-wave im scattering 
lengths [|IJ. This scenario seems to be confirmed by recent Brookhaven data on Kp A decays 
and will be further scrutinized when the pionium lifetime measurements performed at 
CERN H have been analyzed. For the three flavor sector, the situation is, however, less 
clear. Indeed, the observation that m s ~ Aqcd has even led to investigations considering 
the strange quark as heavy (Q,^, |5|. Furthermore, there have been recent speculations that 
the structure of the QCD vacuum changes dramatically with increasing number of flavors, 
e.g. it could be possible that the condensate is sizeably suppressed in SU(3) as compared 
to SU(2) 0. 

One of the cleanest processes to test our understanding of the symmetry breaking pattern 
in the presence of strange quarks is elastic pion-kaon (ttK) scattering near threshold. This 
reaction is interesting for a variety of reasons. First, it is very similar to tttt scattering 
in the two-flavor sector but also different in that the quark mass difference m u — md can 
appear at leading order in the %K scattering amplitude. Second, there exist abundant data 
from inelastic processes which allow one to extract low-energy characteristics of the ttK 
scattering amplitude by means of dispersion theory. The existing determinations of the S- 
wave scattering lengths are, however, plagued by large uncertainties, see e.g. |§. For recent 
work on combining dispersion relations and chiral perturbation theory, see ||. Third, the 



DIRAC collaboration intends to measure the lifetime of ttK atoms at CERN |10[ which 
gives direct access to the isovector S-wave scattering length. To also pin down the isoscalar 
S-wave scattering length, one would have to measure the 2P — 2S level shift, similar to what 



has been done for the pion-nucleon system at PSI [IT]. The precise relation between the 
7rK atom lifetime and the scattering length (the so-called modified Deser formula) can be 
worked out by means of an effective field theory for hadronic bound states (for the pionium 
case see e.g. [O, [R 



To achieve the necessary accuracy in the ttK system, it is mandatory to sharpen the existing 
one-loop predictions for the scattering lengths and range parameters l| 1 ! ] by including 



isospin violation due to strong and electromagnetic effects. This is done here. First, we 
study isospin violation for the S-wave scattering lengths at tree level for all physical chan- 
nels. This allows for a first estimate of such effects and is also interesting to compare directly 
to the ixix case. Then we focus on the one-loop strong and electromagnetic corrections to 
the relevant channel for ttK atoms, n~K + — » tt°K°, including also a complete treatment of 
soft photon radiation, n~K + — * 7i°K 0, y. 

2 Lagrangians 

In this section we discuss the strong and electromagnetic effective Lagrangians underlying 
our calculations. All the pieces have been discussed extensively elsewhere in the literature, 
so we will just present the terms needed for the following. The effective Lagrangian can be 
expanded at low energies according to 

r _ _ r {2) + £( 2 ) + r (4) + r( 4 ) + m 

*-eft *-"str ' '"em ' '-'star ' '"em ~ ■ ■ ■ j \- L J 

where the superscripts (2), (4) refer to the chiral dimension and "str" and "em" denote the 
strong and electromagnetic terms, respectively. Chiral power counting for the strong sector 
attributes the chiral dimension q to all pseudo-Goldstone boson masses (M^, Mk, M v ) as 
well as all momenta involved. This scheme is generalized to include electromagnetic terms 
by counting the electric charge e also as a quantity of order q, which is dictated by the 
requirement that a unique dimension be assigned to the covariant derivative (that includes 
the lowest order photon coupling). Therefore any term of the form e 2j q 2k M l n M^ M™ with, 
e.g., 2j + 2k + I + m + n = 4 is counted as fourth order. 

The lowest-order Lagrangian is given by 

£ (2) = -\f, v f^ -^{d^f 

+ ^(D IM tfD"U + X U ] + X^U) + C(QUQrf) . (2) 

Here, (. . .) denotes the trace in flavor space. F^ u = d fl A u — d u A fl is the usual electromagnetic 
field strength tensor, A refers to the gauge fixing parameter. All calculations were performed 
in the Feynman gauge, A = 1. U — exp(z$/F) collects the (pseudo-)Goldstone boson fields, 
the low-energy constant (LEC) F can be identified with a common meson decay constant in 
the chiral limit. D^U = d^U — i[v^, U] — i{a^, U} defines the covariant derivative acting on 
U in the presence of external vector (v^) and axial vector (a M ) currents. \ includes external 
scalar (s) and pseudoscalar (p) sources, x — 2B(s + ip), where for our purposes, only the 
quark mass term in the source s is of relevance, s = M. + . . . , M. = diag(m„,m rf , m s ). 
The LEC B is linked to the quark condensate at leading order as already discussed in the 
introduction. Q is the quark charge matrix, Q = ediag(2/3, — 1/3, — 1/3). The LEC C 
accompanying the last term in the Lagrangian eq. (y) can be calculated from the leading 
order electromagnetic pion mass difference, M 2 ±—M 2 = 2Ze 2 F 2 (neglecting a tiny strong 



mass difference ~ (m u — m^) 2 ), where we have defined the convenient dimensionless constant 
Z = C/F 4 . Using F = F w = 92.4 MeV, one obtains Z « 0.8. 



The fourth-order Lagrangian [I5|] contains the following "strong" terms which are needed 
for the amplitude describing the process ir~K + — > tt°K° at one-loop level: 



(4) 

+ L 5 (D^D^U^U + tf X )) + U (x ] U + x^f 

+ L 7 (x'U - X U ] ) 2 + L 8 tfU^U + xrfxrf) • (3) 



C% = L 3 {D^D»UD v tfD u U) + L A (D^D^U) { X ] U + X U ] ) 

, I'll' I ! 1 / /.■'»' I — i— .' X ' ' I XI 4 



With regard to the analytic formulae given for the S-wave scattering length in app. |A], 
we remark that only three of the LECs defined in eq. ([|) play a role for this quantity at 
leading order in isospin violation. Terms depending on the Zweig rule suppressed constants 
Z/4 and Lq as well as on L 7 are of higher order in isospin breaking. In addition, the following 



electromagnetic counterterms given in |16| are needed: 

£<£ = K 1 F 2 (D ll U^D^U){Q 2 } + K 2 F 2 (D fl U^D^U}(QUQU^} 

+ K 3 F 2 ((D I1 WQU)(D^WQU) + (D (t UQrf)(D»UQrf)] 

+ K 4 F 2 (D^QU)(D fl UQU^ + K b F 2 {(D ^ D»U + D^UD^)Q 2 ) 

+ K 6 F 2 (D^D^UQWQU + D fl UD"U*QUQrf) + K 7 F 2 ( X ] U + X U ] ) (Q 2 ) 

+ K 8 F 2 ( X ^U + xU^)(QUQU^) + K 9 F 2 ((x ] U + \J\ + xU ] + U X ^Q 2 ) 

+ K w F 2 ((x ] U + U^x)QU ] QU+ (xW + Ux^QUQrf) 

+ K n F 2 ((x ] U - U^x)QU^QU + ( X rf - Ux ] )QUQU^) 

+ K 15 F 4 (QUQU^) 2 + K W F 4 (QUQU^(Q 2 } . (4) 

Also the dependence on a few of these (K 7 /g, Ki5/ie) cancels in the S-wave scattering 
length. One more term (proportional to K12 in ||16|| ) that is needed in principle for the 
renormalization of the decay constants of charged mesons can be omitted in case one uses 
physical values for decay constants with electromagnetic effects already subtracted (see 
section 14.1]) . 



For numerical evaluation, we use the central values and error estimates for the hadronic 
low-energy constants as given in [I7|. For the electromagnetic ones, we use the estimates 
obtained via resonance saturation in fL8 |, and add error bars of natural size (±l/167r 2 ) 
uniformly. 

3 Isospin violation at tree level 

Pion-kaon scattering in the isospin limit can be described by two independent amplitudes 
T 1 / 2 and T 3 / 2 , corresponding to isospin 1/2 and 3/2, respectively. It is sometimes convenient 
to combine these into isospin-even and -odd amplitudes T ± which are defined by 

T a(i = 5 a pT + + - [r Q , Tp}T~ (5) 



(a, (3 refer to the isospin indices of the pions) and which can be related to the amplitudes 
of definite isospin via 

3T+ = T 1/2 + 2T 3/2 , (6) 

3T~ = T 1//2 — T 3//2 . (7) 

All these amplitudes can be expressed in terms of the usual Mandelstam variables s, t, u. 
They can be decomposed into partial waves tj(s) using 

T 7 ( S ,£) = 167r]r(2/ + l)£f( S )P^) , (8) 

l 

where z = cos 6 denotes the scattering angle in the center-of-mass system, and Pi{z) are 
the Legendre polynomials. Close to threshold, one can parameterize the real parts of the 
partial wave amplitudes in terms of scattering lengths (ai) and effective ranges (bi), 

Retf(s) = Y- (|q in ||qout|y {a/ + &f|q in | 2 + O (|q in | 4 )} , (9) 

where we have already accounted for the possibility of different masses of the incoming 
and outgoing particles and therefore for different incoming and outgoing center-of-mass 
momenta, 



v (s - (Mf-Mff)(s - (Mj? + Mr) 2 ) 
lq.nl = ^ , (10) 

Ms - (Mf-Mf) 2 )(s - (M^+¥f) 2 ) 

iQoutl = ^= , (11) 

where M™™ refer to the physical masses of the incoming and outgoing pions and kaons, 
respectively In the isospin limit, eq. @ collapses to the usual definition. We note that 
the real part of a total amplitude at threshold (i.e. for |q in | = 0) is linked to the S-wave 
scattering length by 

ReT/ hr = 87r(il4 n + ilC) a T + O (|q in | 2 , |q in ||q out |) . (12) 

At tree level, the isospin-even and -odd pion-kaon scattering lengths are given by [|19[ ^ 

a+ = , a^ = M " Mx = 70.8 x 10~ 3 /M W . (13) 

° 8ttF2(M^ + M x ) ' K J 

The different physical pion-kaon channels receive corrections to their isospin symmetric 
scattering lengths when taking into account the mass differences as well as insertions stem- 
ming from the electromagnetic term in eq. (0). We collect the isospin breaking parameters 
as S G {m u — md,e 2 } and expand the corrections to the scattering lengths up to order 5. 
The following conventions were used: 



The isospin symmetry limit is defined according to e.g. [1^1 , i.e. we express everything 



in terms of the charged meson masses M„±, M%±. Note that these are not exactly 
the natural choices from the point of view of a chiral analysis of the meson masses, 
in which case one would prefer to take the neutral pion mass as a reference (which is 
indeed done in PH ), as one has M 2 = 2Bm (with rh = ^(rriu + md)) to good accuracy, 
neglecting only tiny corrections of order {m u — m^) 2 . However, one would consequently 
have to resolve to using a "non-physical" kaon mass M\ = B(fh + m s ) ~ 495 MeV. 
Arguments for the use of the charged pion and kaon masses are the correct kinematics 
for the experimentally accessible channels (which use incoming charged particles), 
and the fact that the existing study of pion-kaon scattering in ChPT (in the isospin 
symmetry limit) |14| also employs these choices for numerical evaluation. 

As finally we want to evaluate isospin violating effects up to order e 2 and rriu — md, 
we can neglect the strong pion mass difference (of order {m u — m d ) 2 ). Quark mass 
insertions can then, at leading order, be expressed by physical meson masses according 
to 

M 2 M 2 o - M 2 ± + M 2 ± - M 2 



m — > — — , rrid — m. 



w , B 

M^+M^-M^ M 2 ± -M 2 q 

One important isospin breaking effect due to the light quark mass difference is 7r°r]- 
mixing: the 7r° and rj fields are given in terms of the SU(3) eigenstates cf> 3 , cf) 8 according 
to 

*°) = ( C ° Se Shle )(t 3 ), (15) 

T] J V — sin e cos e / \ 0g / 

where the 7r°?7 mixing angle e is given by 

1 (V3m d -m u \ 

e = - arctan . (lb) 

2 \ 2 m s — m J 

Replacing the quark masses by meson masses according to eq. flj~4]) , one finds the 
numerical value e = 1.00 x 10 -2 . We find it convenient to express all corrections due 
to the light quark mass difference in terms of this mixing angle e, which is of course 
of order m u —md. 

Note finally that we only display scattering lengths for processes involving kaons of positive 
strangeness (K + , K°) as the strong and electromagnetic interactions obey charge conjuga- 
tion invariance, such that scattering lengths for all channels involving scattering of K~ or 
K° can be obtained from those given below. We find: 



3/2 



{l - 0.018} + 0{5 2 ) , (17) 



ao( 



> { ° ° J \ V3M 2 K Ml(M K +M n ) 

= (a+ + oo ) {l + 0.002 - 0.001} + 0(5 2 ) , 
nU+K°-+*»K+-\ fin-U M 2 K -2M K M 7T + 2M 2 Z(M 2 K + M 2 )F 2 

On 7T K — > 7T K = V 2 a n < 1 p= € —p. — 7-— — — , 

+0{5 2 ) 

= V2a {l -0.003- 0.008} + 0(5 2 ) , (19) 

M 2 
a (7r°K + ^n°K + ) = a+ = ^* e + 0{8 2 ) 




c/ 







2.9 x 10" 3 / M - + 0{5 2 ) , (20) 



ao(n K + — > n K~ 




2ZF. 





M^Mk 
= (a+ + an-){l + 0.018} + 0(6 2 ) , (21) 

= -v 7 ^ an"{l + 0.006 + 0.009} + 0(5 2 ) , (22) 

M 2 

= 4 + 2.9 x lO -3 /^ + 0(5 2 ) , (23) 

= ao /2 {l + 0.002 - 0.001} + 0(5 2 ) . (24) 

For the elastic charged channels, eqs. fll7f) , (|2l|) , subtraction of the one-photon exchange 
Born term is required in order to arrive at the above results (see |22| for the comparable 
mi case). These charged-particle Coulomb interactions in principle require a separate 
treatment tailored to scattering or bound-state problems. (For the txtx scattering problem, 
this is discussed e.g. in ||23| .) We remark that isospin violation effects can only be given 
in absolute size for those amplitudes which are (in the isospin limit) proportional to T + , 
as 4 vanishes at tree level. The one-loop corrections to at are rather large, a = 
31.9 x 10 _3 /M 7r |L4]], but still then the corrections displayed in eqs. ( pOf ), (|23|) amount to 
10% effects. (The caveat is, of course, that isospin violating contributions at one-loop level 
might reduce the leading-order terms considerably.) In fact, this is neatly comparable to 
pion-nucleon scattering: as pointed out by Weinberg, the isoscalar ttN amplitude vanishes 
at leading chiral order |H| and is therefore prone to display large isospin violation effects 
24fl , which have indeed been confirmed in [p5fl . 



It is instructive to compare the results above also to the corresponding corrections calculated 
for mi scattering. First of all, m u — md effects are suppressed in SU(2), such that the 



in 



only isospin breaking corrections of interest in the titi case are of electromagnetic origin. In 
contrast, strong isospin breaking appears at leading order in tcK scattering and, as shown 



above, is potentially of equal significance as the electromagnetic effects. On the other 
hand, the authors of |22| find corrections at tree level of the order of 6%, while the largest 



corrections given above are 1.8% (for the elastic charged meson channels). This is not due to 
a particular "smallness" of e.g. the electromagnetic contribution for 7rK scattering, but due 
to the fact that the (isospin conserving) S-wave scattering length is larger, a (iiK) ~ M^Mx 
in contrast to a (7T7r) ~ M 2 . As already hinted at in the introduction, processes with a 
conserved number of strange quarks (such as irK scattering) can also be studied in a chiral 
SU(2) theory in which Green's functions are only expanded around the limit m u = m^ = 0, 
while m s is held fixed at its physical value. Consequently, only the pion is considered to 
be "light" , while the kaon is now a heavy particle and has to be treated in analogy to e.g. 
nucleons in ChPT, resulting in what one might call "heavy-kaon" ChPT [EL El. The ttK 
Lagrangian then starts out at order q, as does the (isovector) irK scattering length (see 
again the analogy to ttN scattering). In such a counting scheme, electromagnetic corrections 
to ttK scattering lengths are therefore necessarily suppressed by one chiral order (as they 
scale like e 2 ), which sheds some light on why these corrections are smaller here than for tttt 
scattering. 

We want to emphasize that it is in general impossible to arrive at the results above by assum- 
ing the dominance of certain "effects" for the violation of isospin symmetry, e.g. studying 
the implications of meson mass splittings or 7i°r] mixing alone. From the point of view of 
ChPT, such a partial analysis is not very useful and can lead to strongly misleading results. 
In order to demonstrate this, we show such a decomposition (artificial though this really is) 
for the specific channel studied in more detail in the following chapter, tt~K + — > tc°K . We 
distinguish "kinematical" effects (due to meson mass splittings, entering essentially via cor- 
rections to the values of u and t at threshold), "7r°?7 mixing" effects which modify the isospin 
symmetric amplitude by factors of sin e or cos e, "quark mass insertions" for the four-meson 
vertex (which vanish in the isospin limit for the channel in question), and "electromagnetic 
(em) insertions" (~ Z) for the four-meson vertex. The relative corrections can then be 
decomposed as follows: 

e e L M K 2M 2 K + Ml M 2 K - Ml 



strong : -== = -^=^1 — ^1^ + K " + -^ -^ , (25) 

6 V3 V3{ M n 3M K M W 3M K M n J ' V ) 

kinematical tt^tj— mixing quark mass 

Ze 2 F 2 Ze 2 F 2 Ze 2 F 2 

electromagnetic : — — -£- = — — - — -^ — h — — - — -^- . (26) 

6 M K M n 2M K M n 2M K M n v ; 

V v ' V v ' 

kinematical em insertions 

It is obvious from the above that for the strong isospin violating contributions, individual 
"effects" are much larger than the total sum. 

4 n~K + — > 7T°K° at one— loop level 

4.1 S— wave scattering length up to fourth order 

In this section, we work out the complete one-loop corrections including isospin breaking 
for the pion-kaon channel that is relevant for the lifetime of irK atoms, n~K + — > tt°K°. At 

8 



this level, virtual photon loops play a role and enter the amplitude both via wave function 
renormalization and electromagnetic corrections to the (lowest order) strong vertex, see 
fig. |l|. As is well known, the photon exchange diagram (a) in fig. [I] induces a pole at 




(a) (b) (c) 

Figure 1: Virtual photon loop diagrams for tt~K + — > tc°K°. Crossed diagrams are not 
shown. Diagram (c) denotes effects entering the scattering amplitude via wave function 
renormalization of the incoming charged mesons. 

threshold, the so-called Coulomb pole, which resides in the loop function G nK ' y (s) discussed 
in detail in app. |B|. This Coulomb pole, however, can be calculated and subtracted from the 
amplitude unambiguously. The threshold expansion eq. (|T2p therefore has to be modified 
for tt~K + — > 7T°K° in the presence of virtual photons in order to sensibly define corrections 
to the S-wave scattering length: 

ReT = ^M K ±M ± af + g/ M +M A aQ + c(| () (2?) 

q. v / 

in | 

where a denotes the S-wave scattering length in its tree level approximation (given ex- 
plicitly up to order 5 in eq. @3)), and a is the S-wave scattering length to be defined here, 
including all corrections. 

A further problem occurs when including virtual photon loops: the amplitude becomes 
infrared divergent. We will give details on how to deal with these infrared divergences in 
the following section. For the moment we only note that these vanish at threshold and 
therefore do not affect the definition of the S-wave scattering length as given in eq. fl2?|). 

Analytic formulae for the S-wave scattering length at one-loop accuracy are given in app. |A]. 
In order to arrive at these, the following effects were taken into account, as well as the 
following conventions were used in addition to those already mentioned in the previous 
section: 

• In order to have a strong check on the scattering amplitude, we have checked the 
cancellation of divergences in the exact expressions, i.e. without any expansion in e 
and e 2 . For this purpose, we have used the /3-functions for the various counterterms 
as given in [15, [16 . 



For quark mass insertions at tree level replaced by meson masses according to eq. (|14D, 
the renormalization of these meson masses has to be taken into account properly. 



The mass of the rj is always taken as given by the Gell-Mann-Okubo relation, which 
reads 



?>Mi 



'M 



. „ 2M Z K± + 2M Z K „ - 2M< ± + M^ 

when taking into account isospin violation up to order e 2 and m u — md. Deviations 
from this relation are beyond the accuracy considered in this paper, as the ?7-mass 
only enters via loop diagrams. 





"^0» 



Figure 2: Diagrams contributing to 7T 77 mixing at next-to-leading order. The labels "7r 



and "77" refer to the tree-level mass eigenstates. 
insertion (strong or electromagnetic). 



The black square denotes a fourth-order 



Details on how 
It is probably easiest to 



7T°7/- m ixing j^g t be taken into account at next-to-leading order, 
this is to be done can be found e.g. in |26], |7|> HI IT 
calculate the diagrams in fig. |2| directly in order to account for mixing beyond leading 
order. We have done the calculations using tree-level mass eigenfields for ir° and 
T) throughout. The off-diagonal element of the self-energy matrix (in the basis of 
tree-level eigenfields) may be written as S 7r o T? (g 2 ) = Z^o^q 2 + Y 7V o v , such that Z 7T o v = 
dY* n o v /dq 2 can be thought of as the off-diagonal analogy to the field renormalization 
factors. Denoting the tree-level amplitude for tt~K + — > rjK by T v , we can write the 
additional mixing diagrams as 



T v x M2 J_ M2 x £.0, (q 2 = M 2 



K>-m$ 



S^(g 2 = i(M 2 o + M 2 )) 



(29) 



c (4) 



e- 4 ' thus defined is the correction to the tree level 7r°T] mixing angle up to order 0(q 4 
as given e.g. in 
is finite. 



29| . In contrast to the wave function renormalization factor Z^, e^ 



When expressing the low-energy constant F in terms of physically observable decay 
constants, we always renormalize it such as to obtain F n (as opposed to Fk, say). 
On tree level, the difference between using a normalization of the amplitude equal to 
1/F 2 and 1/F^Fk amounts to nothing more but shifting certain contributions from 
tree to one-loop level, which makes no difference in the sum. However, it turns out 



10 



that the one-loop corrections e.g. for the scattering lengths are much smaller when 
normalizing the tree level expression to 1/F%, which means that in the seemingly 
more "symmetric" case of using 1/F^Fk, the bulk of the (large) corrections at 0(q A ) 
have nothing to do with corrections to the scattering dynamics, but only with the 
renormalization of Fk- 

For the one-loop contributions, the difference between F^ and Fk is of higher order, 
whereas the numerical difference is potentially significant (noting that these scale as 



1/F ). Regarding F ~ 88 MeV |30 as the "natural" constant when writing down a 



scattering amplitude, F w seems more appropriate than Fk- 

Finally heavy-kaon ChPT, which we alluded to before already, also provides some 
insight on F n versus Fk' there is of course no kaon decay constant in a kaon-number 
conserving theory, hence F„ is the only meson decay constant available. As heavy- 
kaon ChPT is supposed to have a better convergence behavior than SU(3) ChPT, 
but ultimately is equivalent up to matching of the different sets of LECs, it seems 
reasonable also from this perspective to only use F v in the nK amplitudes. 

We use the charged pion decay constant in the absence of electromagnetism, as the 
value commonly used for the physical pion decay constant, F n = 92.4 MeV |H|], was 



extracted from charged pion decays with electromagnetic corrections already taken 
care of. The relation to the bare decay constant F is therefore as given in ||15| , and 
independent of photon loop effects or electromagnetic counterterms. 



We display the isospin conserving and violating contributions, at tree and one-loop level, 

( c ?\ 

relative to the isospin symmetric tree level result (which we just denote by a here for 
reasons of brevity). Numerically, these add up as follows: 

a^ = -100.1 x 10^/M"^ , 
a = af ] { 1 + 0.006 + 0.009 

0(e) 0{e 2 ) 

+ (0.121 ± 0.009) + (0.008 ± 0.002) - (0.009 ± 0.008) } . (30) 

V v ' V v ' V v ' ' 

C(p 4 ) 0{p 2 e) 0(p 2 e 2 ) 

All errors quoted for the different contributions above are due to uncertainties in the re- 
spective strong and electromagnetic LECs as given in JT7], [T3|]. The following comments are 
in order: 

• The isospin symmetric corrections at one-loop order are moderate, given the expected 
slow convergence of chiral SU(3). As detailed above, this depends heavily on our choice 
to normalize the tree amplitude by 1/F% and not 1/F^Fk- In fact, expanding the 
one-loop contributions in powers of M n (equivalently to the heavy-kaon approach), 
one finds that they are suppressed essentially by a factor of M% (and not by M n M K 
or M K ) compared to the tree level values, which is what one would expect from an 
SU(2) expansion and which accounts for the smallness of the one-loop corrections. 
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Although they are very small corrections to the scattering length, one might worry 
about the fact that tree and one-loop level isospin breaking corrections (strong and 
electromagnetic effects taken separately) are of equal size. In order to understand 
the relative largeness of these effects, it is instructive to expand the terms given in 
app. [A] in powers of M n . Isospin violating terms are suppressed in comparison to 
isospin conserving ones at the same chiral order by factors which happen to be small, 
but are of chiral order 1. On tree level, these suppression factors for tt~K + — > ir°K 
are e and e 2 /M n Mx, respectively, see eq. fl2"2|). On one-loop level, however, where 
the isospin symmetric loop corrections are particularly small as discussed above, the 
isospin violating ones are suppressed with respect to these only by factors of e M K /M n 
and e 2 M K /Ml (at leading order in M„). 

On the other hand, combining strong and electromagnetic isospin breaking effects, the 
corrections at one-loop level cancel to a large extent, leaving mainly an uncertainty 
due to a lack of sufficiently precise knowledge of the relevant counterterms. This 
cancellation does not take place at tree level. We have checked numerically that no 
significant errors arise from truncating the "exact" expressions at leading orders in 
the isospin breaking parameters e, e 2 . In fact, the corrections arising thereof are at 
least one order of magnitude smaller than the accuracy displayed in eq. 



• With the large number of electromagnetic counterterms, one might have suspected 
that one cannot make any sensible prediction about electromagnetic effects at all. 
While the uncertainty quoted above is as large as the absolute size of the electro- 
magnetic correction at one-loop level, it is, on the other hand, still not larger than 
the uncertainty stemming from our lack of knowledge of the values for the hadronic 
low-energy constants. 

4.2 Infrared divergences and tt~K + — > ir K°j 

As mentioned in the previous section, photon loop corrections induce infrared divergences 
in the scattering amplitude. It is well known how to handle this problem, namely by 
introducing a small photon mass m 7 which brings the infrared divergences into a form 
~ logm 7 . Such terms enter both via wave function renormalization of charged mesons 
and via the loop function G nK ' y (s), as discussed in app. |B|. One then has to include the 
corresponding radiative process (tt~K + — > 7r°K 0/ -f for the case in question) up to some 
maximal photon energy, given either by the maximal energy available due to kinematics, 
or by some experimental detector resolution AE below which soft photon radiation cannot 
be discriminated from the non-radiative process. We are interested in the threshold region, 
where the kinematical limit for the photon energy, 

s-(M K o + M w of 
E, = ^= , (31) 

numerically amounts to 0.6 MeV, such that it seems reasonable to integrate the full cross 
section without making additional cuts on the photon energy. The infrared divergences 
cancel upon adding up the cross sections for both non-radiative and radiative processes. 



This is demonstrated in some detail in app. |C.l . 
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We have emphasized in the previous chapter that the infrared divergences induced by photon 
loops vanish at threshold, such that one might even omit all these considerations when one 
is exclusively interested in the S-wave scattering length. The authors of [^ have argued, 
however, that one should rather define the S-wave scattering length from an infrared-finite 
quantity, which would be the combined total cross section. For an amplitude which can be 
expanded according to eq. (0), the total cross section at threshold is given in terms of the 
S-wave scattering length as 

athr = ^47r(a ) 2 . (32) 



One may now replace the cross section on the left-hand side of eq. (|32[ ) by the combined 
total cross section o + cr 7 and define do via eq. (|32|). 

This is what was suggested in p2| for the similar case of the scattering process tt + tt~ — > 



7r 7r (7), where the authors claim to find a shift in the (redefined) scattering length induced 
by the (infrared-finite) remainder of the radiative cross section at threshold. We have 
however recalculated the cross section for 7r + 7r~ — > 7r°7r°7 at threshold and find that it 
vanishes. This is in our opinion rather obvious for the isospin symmetric case, where the 
neutral pions in the final state are of the same mass as the charged incoming ones. One 
would therefore expect the corrections to the scattering length due to soft photon radiation 
to be small compared to other isospin breaking effects, as they should be suppressed by at 
least two powers in isospin breaking parameters (a factor e 2 in the cross section due to the 
photon coupling, and at least one power of the pion mass difference). The fact that this 
cross section vanishes exactly at threshold in the nir case is due to the enhanced symmetry 
of this reaction, see app. |(J.2| . 

We do not find the cross section for tt~K + — > 7r°7^°7 to vanish exactly at threshold, but to 
be highly suppressed: the leading contribution is of order 6 4 in isospin breaking parameters, 
which is clearly way beyond the level of accuracy we achieve in our calculation of isospin 
breaking effects in the non-radiative amplitude. Details for the calculation can be found in 



app. |C.2| , the result found there is 



In addition to the high power in 5, there is even an unnatural suppression in the mass 
difference (M K ±— M K o+M n ±— M n o) which is due to a numerical cancellation of strong and 
electromagnetic isospin breaking effects. If one calculates a shift in the scattering length 
defined via eq. (^) which is due to the inclusion of the radiative cross section, clq — > a + ag, 
the quantity Oq is found to be 

7 e 2 (M K ±-M K o + M n ±-M n o) 3 5 

00 840V2tt 3 F2 (M K± +M n± y + U[ ° } 

= -2.2 x 10- U /M n ± = 2.2 x KT 13 x af , (34) 

which is a suppression of about eleven orders of magnitude compared to the shifts in the 
scattering length calculated before. We therefore conclude that bremsstrahlung corrections 
are totally negligible for any analysis of tctt or tcK bound state experiments. 
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4.3 S— wave effective range, P— wave scattering length 

Although a precise knowledge of the S-wave scattering length is of highest interest for 7rK 
atom studies, one can of course also study isospin breaking corrections to other threshold 
parameters, particularly to the S-wave effective range b and the P-wave scattering length 
a±. We will not show results for these in as much detail as for a^, but concentrate on the 
channel for which we have calculated the one-loop corrections to isospin breaking. We 
remark that, as indicated in eq. (p7[), photon loops also induce a term linear in |qi n | in the 
threshold expansion of the scattering amplitude, which we disregard here (as it can also be 
calculated and subtracted unambigously) . 

The isospin symmetric tree level results for these two quantities are given by 

b® = -V2b- = Ml + Ml (35) 

° 8x/2vr FIM^M K (M K +M^) ' 



,(2) 



-V2a7 = = . (36) 



The isospin violating effects at tree level can then be expressed as 

*(»-*♦ -. M) = -V2K {! + -£=- j^} + 0(0 , (37) 

ai(jt-K + -^ir°K°) = -V2a^ {l - v^e} + 0{5 2 ) . (38) 

No electromagnetic effects affect the P-wave scattering length at this order. 

We evaluate the fourth-order contributions to these threshold parameters numerically only, 
therefore we do not differentiate between strong and electromagnetic isospin violation at 
this level. Again, all corrections are given relative to the isospin symmetric tree-level result: 

b™ = -54.0 x 10~ 3 /M3 ± , 
b = b { o ] { 1 + 0^06 - O002 + (0.012 ± 0.038) - (0.013 ± 0.005) } , (39) 

af = -4.7 x 10- 3 /M 3 ± , 
a x = a{ 2) { 1 - 0.02 + (0.53 ± 0.14) - (0.03 ± 0.01) } . (40) 

*• "^tC " * ' v * 

°( e ) C(p 4 ) 0(p 2 5) 

We note that the isospin symmetric one-loop corrections to the P-wave scattering length 
are rather large compared to the leading term. The isospin breaking corrections for both 
quantities are of similar size as for ao, only slightly enhanced for ax- The main difference is 
that the cancellation between electromagnetic and m u —md effects at fourth order observed 
for ao does not occur here, such that isospin breaking at fourth chiral order is equally 
important as on tree level. If we compare the one-loop contributions only, the isospin 
breaking effects on the S-wave effective range appear to be very large, equally large even as 
the isospin symmetric one-loop corrections. However, as already hinted at by the error range 
for the latter, this does not result from an unnatural enhancement of isospin violation effects, 
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but from an accidental suppression of the isospin conserving part. Individual contributions 
to the latter (from separate loop diagrams, counterterms) are much larger. 

In contrast to both the S-wave and the P-wave scattering lengths, the S-wave effective 
range as defined naively from the scattering amplitude alone is infrared divergent. The 
above results were achieved by setting m 7 equal to the maximal energy of a bremsstrahlung 
photon, which at threshold is roughly 0.6 MeV. This is what effectively happens when ren- 
dering infrared divergent expressions finite by including the appropriate radiative process, 
while we neglect here any finite contributions from the latter, see the discussion in the pre- 
vious section. We also note that a redefinition of these higher-order threshold parameters 
from the infrared-finite total cross section (as done in the previous section for the S-wave 
scattering length) cannot be done unambiguously. 

5 Summary 

In this work, we have considered isospin violation in pion-kaon scattering near threshold. 
To systematically account for such effects due to the light quark mass difference as well as 
electromagnetic interactions, we have made use of SU(3) chiral perturbation theory in the 
presence of virtual photons. The pertinent results of this investigation can be summarized 
as follows: 

(1) Already at tree level, one has strong as well as electromagnetic isospin violation. We 
have considered all physical channels and found that these effects are in general small 
(at most two percent). In particular, the relative corrections are smaller than in 
the comparable case of mi scattering where isospin violation at tree level is a purely 
electromagnetic effect. We have also stressed that considering one particular source 
of isospin violation only can ensue very misleading results. 

(2) Because of its relevance to the lifetime of ttK atoms, we have considered the one- 
loop corrections to the S-wave scattering length for the process n~K + — »■ tt°K°. 
The fourth-order isospin symmetric corrections are moderate and the isospin vio- 
lation effects at this order cancel to a large extent (for the central values of the 
low-energy constants used here). The uncertainty in these isospin violating contri- 
butions is mainly due to the poor knowledge of the electromagnetic LECs K t . This 
uncertainty is comparable to the one in the isospin symmetric amplitude induced by 
the variations in the strong LECs Lj. 

(3) We have also considered the radiative process tt~K + — ► 7c°K 0/ -f. We have explicitly 
demonstrated the cancellation of the infrared divergences. The remaining radiative 
cross section is strongly suppressed at threshold so that the properly redefined (finite) 
scattering length acquires no significant correction. 

(4) We have also given the one-loop corrections for the effective range bo and the P-wave 
scattering length Oi for this channel. The isospin violating corrections are somewhat 
more pronounced than for the S-wave scattering length. 
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As pointed out in the introduction, for precisely predicting the 2P — 2S level shift in the 
■kK atom, one would have to extend these considerations to the elastic scattering channel 
n~K + — > 7i~K + beyond the tree level result presented here. 

After submission of this work, a similar analysis with comparable results has appeared, 



see |2] 



We thank G. Ecker, A. Nehme, and P. Talavera for helping to resolve an error in the first 
version of the manuscript. 

A S— wave scattering length at fourth order 

In this appendix, we show the analytic expressions for the one-loop corrections to the scat- 
tering length ao for tt~K + —>■ n K°. The isospin violating effects are expanded up to order 
e and e 2 . We express all these as relative corrections A^Oo/cio • The low-energy constants 
with the infinite part subtracted are denoted by L\, K?, as done conventionally. Note how- 
ever that we do not display the dependence of these various constants on the renormalization 
scale A explicitly. All corrections given below are of course scale independent. 

The isospin symmetric contribution is 
AgU _ MM 



c*6 

1 ( 1 



16tt 2 \M 2 K -Ml 



4M K 
+ 



4MJ - |M|) log f - f Ml log f + (£4 - f ) log 



9M OT 



(M K - M„){2M K + M w ) f2(M K + M 77 ) I M K - M w 

arctan - 



M K + M n \ M K - 2M n V 2M X + M„ 



(M K + M„)(2M K - M v ) f2(M K -M n ) I M K + M/ 

M K -M V \ M K + 2M n V 2M K - M, # 

The strong isospin violating part can be expressed as 
A 2°o e f^,,,,, SifW , 8(2Mi-M|M 2 -3M i ,M w 3 -M^ 



(41) 



16M^(M i ,-M^)L 3 + 



af V3F* { * v A " 7 ^ 3M^ 

16(Mj-M 2 )(2M| r + M 2 ) r 



1 f72Mi - 56MiilC - 



167T 2 



72M 5 K - 5QM^M n - 20OM|jg + 2\M 2 K Ml - 19M X M^ + 2Ml M^ 

QM K (M K - M n ) 2 (M K + il4) * ° g A 

94M£ - 185MIM,- - 915Mg-M 2 + 46M K M* - 20M* M K 

27(M K - M v y(M K + M v ) * ° g A 

/ A8M 7 K - A%M & K M n + 388M|-M2 - 224M£il^ 
V 54M K M W (M K - M„)*(M K + M„) 
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+ 



192M&M* + 2h9M 2 K Ml - 269M K M% - 6M 
+ 54M K M n (M K - M V ) 2 (M K + M, 

2M 7r (16Mg r - 40M K M W - 47M 2 ) 
9(M K - M n ) 




2M K M n (4M K + M n ) / Ma - Mr / 2(Ma + M„) M K - M w 

3(M K + M w ) V 2Ma + M n arC an ^ Ma- - 2Af w V 2M K + Mr 
2M^(116M^ - 79M\M^ - \5QM 2 K Ml + 74M K M^ + 8M*) 
27(MA-M^)2(MA + M 7r ) 



M K + M n 2{M K -M V ) M K + M„ 

V 2M X -M^ V M K + 2M n V 2M X - M ff# 

Finally, the electromagnetic corrections at one-loop level can be written as 
AgS = 2\ Z ( SL *(2M K -M„)(2M K + 7M n ) \ 



(42) 



*W j_ *-n -l. 6M ^ + 2Mj^ - 6MaM 2 + Ag 

Ma(4M| - 7M, 2 ) 2 hnjr+jm 

9M.(M|-M 2 ) (^5+ ^)" 9(10^5 + ^) 

_ M 7r (2Mj + M 2 ) 2(4Mj - M 2 ) \ 

9M K (Ml-M*) [K9+Kw)+ 3M K M n ^ K io + K n)j 

Ze 2 f 40M| - 118M^M^ + 22M|-M 2 + 26M 2 r M 3 - 71M K M* + HMg M, 
+ 16tt 2 \ QM K (M K - M K )\M K + Mr) ° g A 

36M|- - 97M^M W - 211M|-M 2 + 126M|-M3 - 400MaM^ + 56M| Ma 
27M w (Ma - M 7r ) 3 (MA + M T ) ° g ^T 

32M| + 296M^M 7r - 270MIM 2 + 2^0M|M| - 121MaM^ - 27M^ M, 
54Ma(Ma - M n ) 3 (M K + M^) ° g T~ 

113M£ - 138M|-M^ - 249M X M 2 - 13MaM 3 + 3M^ 
18M K (M K - M n ) 2 (M K + M n ) 
16M 2 - - 17M K M n - 8M 2 f2(M K + M n ) I M K - M n 



+ 



arctan 



9(Ma + M w ) ^{2M K + M^) (Ma - M ff ) \ M K - 2M n V 2M X + M n 

32M K + 35M|M^ - 108MIM 2 - 40MaM 3 + 44M^ 



27(Ma - M 7r ) 3 A /(2MA - Mr) (Ma + Mr) 



+ {2(M K -M„) Ma + M. 
x arctan 




2M K -M« 



e 2 f 12M 2 - - 4M X M, - M 2 M„ 4M 2 - + MaM, + 12M 2 , Ma 

log — - — I — — — - log — — 



16tt 2 [ 2Ma(Ma + M^) & A 2M n (M K + M n 

4M 2 - - 18M K M W - M 2 
3M K M n 
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(43) 



B Infrared divergent loop diagrams 

The essential loop function containing an infrared divergence, stemming from the photon 
exchange diagram (a) in fig. [I], is 



G KK ' r {s) = -i f 



d A k 



( 27 4 ((q, + k) 2 - M 2 ± ) ((q K -k) 2 - M 2 K± ) (k 2 - m 2 ) ' 



(44) 



where s = (qx + Qv) 2 - Evaluating this in the kinematical region needed for the process in 
question, s > (M K ± + M n ±) 2 , we find 



Re G wK ~<(s) 



32-K 2 as 



log a — log 



m; 



1 1 — ^1 , i Z 2 

log h log 



+- 



Z\ 



l-z 2 



(45) 



log (l—zi) - log z x + log z 2 - log (l-z 2 



i i Z2 \ Zl 
+ log — log — 

Z 2 



a 



+ log log 1 



a 



1-2, 



^2 



Li 



(7 



l-*i 



4tt z 



(9(m 7 ) 



7/ ' 



Im G^ 7 (s) 
where 



-J_{bg('^)-21og«T^ + C)(m., . 

l07T<7S L \ S 



(46) 



1/ A\ a 

Z ^ = 2{ l --) T 2' ° 

S = M 2 K± +M 2 ± , A = iU 



2S A 2 



and 



Li(z) 



log* 



dt 



(47) 



(4£ 



'i 1-t 

is the dilogarithm or Spence function. We have checked that this loop function coincides 
with the analogous one needed for 7T7r-scattering as given in [[J3|, [22] when going to the limit 
of equal pion and kaon masses. Other infrared divergent loop contributions stem from the 
wave function renormalization of the charged pions (see e.g. ||34|| ) and kaons, 



Z K ± 



,^ log ^L + OK) 



7 y > 



i lo ^ + « 



We write the tree level amplitude for n K + — > 7i°K° as 

T (2) = c + c s s + c t t + c u u , 

where 

B ( sine\ Ze 2 ( sine\ 

c ° = a /2F2 1 ( m u- m d) cose + {m u + m d -2m s )—^j y=\ cose + —=■ J 



2y/2F 2 



cose + 



sine 



y/5J' 



Ct 



sine 

Vqf 2 



1 



2V2F 2 



cose 



sine 
75 



(49) 



(50) 



(51) 
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The infrared divergent piece of the real part of the one-loop amplitude T^ is then given 

by 



(ReT< 4 > 



div ,„,. e 2 f s — S 

= T (2) x — 



7r" I as 



1 1_21 , 1 Z2 
log h log 

z\ l-z 2 



= T(2) x^K)- 



log 7i~ log ii4r" log M^I 

(52) 



The total cross section is obtained from the amplitude via 

1 |q ou t| 



a 



647r 2 s Iqj, 



\T\ 2 dtt 



(53) 



with |q in | and |q ou t| defined as in eqs. (PUD, (|TT|). The infrared divergent part of the total 
cross section can be calculated from eq. (|53|) to be 



a 



div 



Co + c 



Q + C u \ Cj ~\~ C u / 2 



s + 



(M| ± +M|-o + M 2 ± +M 2 o) 



e I qout I 

(4vr) 3 |q in |s 

c t -c u (M 2 :± -M 2 ± )(M| r o-M 2 o)l 2 j _4 2 2 2 \ 

+ -(Q-c u ) |q in | |q out | ^L(m 7 ) . (54) 



C The radiative cross section 




(a) 



(b) 



(c) 



Figure 3: Diagrams contributing to ix K + — > n K 7 at tree level. 

The amplitude for the process 

ir-(q^)K + {q K ) -► TT ^) K°(p K )j(l) (55) 

is given, to lowest order, by 

ee*-{2q K -l) ( >. 

i ' = -7; : 5- c + C s S 2 + Q U + c u u n ) 

2q K ■ I — m* v y 

e e* ■ (2<2V - /) / \ 

- -z j~2 — — \ c o + c s s 2 + c t t K + c u u K J 

+ 7= 2 e* • I (?* - <h + Pk ~ Pn) cose + (q K - q w - 3p K + 3p^j —= 1 , (56) 



2y/2F 2 



V^i 
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where 

S 2 = {Pn + Pk) 2 , U/K = (q-K/K ~ Ptt/k) 2 , U^/K = {q^/K ~ PK/n) 2 ■ (57) 

The three terms in eq. (|56|) correspond to the three diagrams in fig. [j| The total cross 
section can be calculated from this amplitude by 

o 1 = , \. 71 1 , r- I dE i I d %, I dtt Pwl —, ^W l \ rrl^ 7 ! 2 , (58) 

(2^4\ctn\V^J J P J ^8(^o + Ko(l + ^Lcos^))' ' l ' 

where 8 P7r refers to the the angle between the outgoing pion and the axis of the incoming 
particles (in their center-of-mass system), and 9 Plv i denotes the relative angle between the 
outgoing pion and the photon. 

C.l Infrared divergence 

One important check resulting from the exact calculation of the radiative cross section is 
that the infrared divergences indeed cancel when combining radiative and non-radiative 
cross sections, we therefore demonstrate how to isolate the infrared divergent parts in the 
cross section eq. (|58|). These stem exclusively from the lower integration bound of the E\- 
integration, therefore it suffices to collect the leading powers in \jE\ of the integrand. The 
squared matrix element simplifies considerably in this approximation, one may set Si — > s, 
tn/K -^ t, u w /k — > u, and obtains 

I (I ■ q„) 2 (I ■ q K ) 2 l-qn 1-qK J v J 

As is well known from Quantum Electrodynamics, infrared divergences arise only from 
diagrams with soft photon radiation from external legs (diagrams (a), (b) in fig. |3|), and 
indeed contributions from the diagram where the photon is radiated from the 7ri^-vertex 
(diagram (c) in fig. |3|) play no role here. Furthermore, we can set |p„-| = |q ou t| (which is, 
as above, the modulus of the outgoing momentum of the non-radiative process), and use 
the abbreviations y = cos 6^/, z = cos 6^, to obtain 

e 2 la J r E T* r 1 r 1 r 2lT dd , 

a 7 = e jqput /■ -< dEiEi dz dy I a ^i 



(47r) 3 4|q in |s Jm 7 J-i J-i Jo 2tt 

( M 2 ± M 2 ± s-M 2 ± -M 2 K± } / n2 _ / nx , s 

Note that in this approximation, the term Co + c s s + c t t + c u u depends, of all integration 
variables, only on z (via t and u). The Pvr ; integration is relatively straightforward when 
using 

I ■ q n/K = Ei |q in | (SV/a- =F cos a) , (61) 

where 

s T M%± ± Ml ± r - r 

Sn/ K = — n , cosa = yz- Vl-yWl- z 2 cos4> P7r i , (62) 

^IqinlV S 
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such that, for example, the first term in the curly brackets in eq. ( |60|) yields 

M 2 f 1 . ^ dfa 1 M 2 ± yi j^-yg _2_ 

"* -/-i y -/o 2tt (/ • g.) 2 " ^|q in p U \si - 1 - 2S„yz + y 2 + z 2 ) 3/2 " Ef ' 

(63) 
and the other terms can be evaluated similarly The only remaining z-dependence is the 
one which is also inherent in the integration of the non-radiative cross section, 

/ dz(co + c s s + c t t + c u uj . (64) 

Altogether, one ends up with 

7 __ e 2 |q out | {E%**} \s-M 2 K± - M 2 ± (,_ S K + 1 , ,__ S v + 1 







Ct + Cu 



(M 2 K± +M 2 K0 + Ml ± + Ml») (65) 



^ - c u (M 2 K± - M 2 ± ) (M 2 K0 - M 2 



+ g(ct - c u ) 2 |q in | 2 |q out | 2 1 + 0(m°) , 



2 s 

which cancels the m 7 -divergence in eq. 

C.2 The radiative cross section at threshold 

The other aspect of the radiative cross section on which one would like to have analytical 
information is the remainder at threshold. As the incoming charged particles are at rest 
at threshold, the angular integration for the total cross section simplifies considerably: the 
integration dQ Pw is trivial, and so is the integration d<p P7r i, hence the only angular variable 
to integrate is dy = d cos 9 Pv i, which is most conveniently done in the center-of-mass frame 
of the outgoing tt°K° system. Large cancellations take place, such that one obtains a 
surprisingly simple expression: 

e 2 ( R . ^(M K ± + M n ±) 2 

• i sine 



/ dy |T°'|? VnT . = — — jfcose — v3siL~ , 

7-i y| lthr 24F 4 ^ ) M K ±M n ± 

[s 2 - (M K ±-M^) 2 )(s 2 - (M K o-M n o) 2 )(s 2 - (M K o + M n o) 



< 2 



(66) 



In order to obtain eq. (B3), it is essential to take the radiation from the four-meson vertex 
(diagram (c) in fig. [3]) into account, the vertex of which is linked to the n~K + —> 7i°K° 
vertex by gauge symmetry. In fact, we find that just by expressing the coefficients for the 
5-point vertex (the third term in eq. ([56])) by c t and c u , any dependence on cq and c s drops 
out completely. This is why for the case of 7r + 7r~ — > 7r°7r°7, the cross section vanishes at 
threshold: the lowest-order 7r + 7r~ — > tt°tt° vertex is proportional to s— M 2 (in the a gauge), 
and the 7r + 7r _ — > 7r°7r°7 point vertex vanishes (again in the a gauge). Furthermore, the 
result of eq. ( |66|) is proportional to (c t — c u ) 2 which demonstrates that only the particular 
"asymmetry" of the reaction in question allows a finite remainder to survive at threshold. 
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Transforming the remaining integration over the energy of the outgoing photon into an 
integration of s 2 , we obtain the total cross section at threshold via 



7 _ " I "lout I / , 

"*" - (8tt)3 |q in | s L 



1 _|qoJ [' d ( s ~ s 2 \ (fg - (M K 0-M w o)*)(s 2 - (M^O+M^o) 2 ) 
3 |q in |s V S2 { s 2 )\ (s-(M^o-M,o) 2 )(s-(M^o + M,o) 2 ) 

j\\T J \L , (67) 



where S2 hr = {M K o + M n o) 2 is the minimal energy (squared) for the final tt°K° system. 
Eq. (|67D is evaluated to leading order in isospin violation parameters, yielding the result 



In terms of fundamental isospin breaking parameters, the meson mass difference in eq. (|68 
can be expressed as 

B(m u -m d ) 2 2 / 1 1 , 2 



Mk*-M k „ + M,*-M«, = -J Mk " + Ze>F^— + -) + 0(6>). (69) 
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